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Abstract 



A brief overview is given of recent developments and fresh ideas at the in- 
tersection of VT— and/or CVT— symmetric quantum mechanics with super- 
symmetric quantum mechanics (SUSY QM). Within the framework of the 
resulting supersymmetric version of CVT— symmetric quantum mechanics 
we study the consequences of the assumption that the "charge" operator C is 
represented in a differential-operator form of the second or higher order. Be- 
sides the freedom allowed by the Hermiticity constraint for the operator CV, 
encouraging results are obtained in the second-order case. In particular, the 
integrability of intertwining relations proves to match the closure of our non- 
linear (viz., polynomial) SUSY algebra. In a particular illustration, our form 
of CVT — symmetric SUSY QM leads to a new class of non-Hermitian poly- 
nomial oscillators with real spectrum which turn out to be VT — asymmetric. 
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1 Introduction 



The recent growth of interest in the possibility of working with non-Hermitian 
observables in quantum theory (cf. the concise review papers [1]) is mainly 
due to the influential Bender's and Boettcher's letter [2] where its authors 
observed that the spectrum of certain Hamiltonians H ^ H* seems real and 
discrete and bounded below. 

They conjectured that such an observation may find a firmer mathe- 
matical background and explanation in the symmetry of their models with 
respect to the combined action of the parity V and time-reversal (i.e., com- 
plex conjugation) T. This inspiring idea has been further developed and 
re-formulated as proposals of the so called "PT— symmetric quantum me- 
chanics [3], pseudo-Hermitian quantum mechanics [4] and CVT — symmetric 
quantum mechanics [5] . They all deal with more or less the same class of the 
specific non-Hermitian models characterized, in the language of the latter 
reference, by another symmetry operator C which is very conveniently called 
" charge" . 

There exists an extensive literature on PT-symmetric quantum mechan- 
ics [6, 7]. In particular, in a number of papers [8, 9, 10, 11], unexpected 
consequences of the non-Hermiticity of Hamiltonians have been noticed to 
emerge after its supersymmetrization a la Witten [12]. In terms of local 
models H = p 2 + V(x) on the real line (x G 1R) where V(x) = V*(—x), these 
Hamiltonians satisfy the intertwining relation 

H ] V = VH. (1) 

Such PT-symmetric Hamiltonians may have either complex, or real spectra. 
When the VT symmetry remains spontaneously unbroken and all the spec- 
trum is real [2], one has elaborated the concept of quasi- Hermiticity of the 
Hamiltonian [13, 14]. This means that the intertwining relation (1) holds 
also with V replaced by a positive-definite operator = 0^ > which plays 
the role of a metric operator. The physical interpretation of such models 
is standard [15]. When the spectrum is complex, relation (1) can still be 
written with V replaced by a pseudo-metric [16, 17, 18]. 

We shall now generalize the previous considerations to a new type of 
symmetry. The framework of our constructions proposed in our recent let- 
ter [19] will incorporate Hamiltonians with both real and complex spectra. 
Correspondingly, we shall also deal, in general, with non-positive metric (i.e., 
pseudo-metric). As for the case of VT symmetry, the interpretation of this 
type of quantum mechanics can be disputable [14] and might require some 
innovation. However, we stress that, in our framework, we can find models 
which have real spectra, where, in particular, a non-Hermitian Hamiltonian 
is related by similarity transformations not only to a Hermitian operator, 
but, more specifically, to a Hermitian Schrodinger operator. Thus, for these 



2 



cases, we recover the standard quantum mechanics, after similarity. There- 
fore, these cases are not disputable in their interpretation. From a conser- 
vative point of view, one might restrict the interest of our supersymmetric 
approach insofar as one takes it instrumentally as a strategy to find com- 
plex Hamiltonians with real spectrum that do not satisfy VT invariance (see 
Section 4 below). 

1.1 SUSY intertwining relations 

In the same spirit as in Ref. [19], we shall study the intertwining relations 

TE^ = ET (2) 

mediated by the Hermitian operator 

T = CV {T= , (3) 

where V is the parity operator, and C a generalized "charge" operator, as- 
sumed to be a polynomial in the differential operator d/dx. For any Hamilto- 
nian E, Eq. (2) is equivalent to CVT conservation, with T the time reversal 
operator, 

CVTE = ECVT . (4) 

In this paper we shall not discuss in detail the metric interpretation for J 7 , 
but only stress the fact that, if T and E satisfy Eq. (2), then also JF^ 1 (if it 
exists) and E meet an intertwining 

tftjT-i = jr-ijj ; ( 5 ) 

which means that E is pseudo-Hermitian with respect to JF -1 , i.e., 
T~ x — pseudo-Hermit ian [4]. 

This observation may be useful for implementing the metric based on JF -1 , 
when JF" 1 has a better behavior than JF, e.g. with respect to boundedness. 
Nevertheless, in our text we also use for a T satisfying Eqs. (3), (2), the 
word "metric" operator. In fact, Eq. (2) implies that ET is Hermitian. As 
a consequence of (5), if | 0) and | ip) are two arbitrary vectors of the Hilbert 
space L 2 (R), we have 

J <t>*{x) (F^Eip) (x)dx = J ip*(x) (F-*H<f>) (x)dx. 

This can be interpreted as a Hermiticity condition for E provided the scalar 
product is defined as 

(0 | = / <t>*{x) (^V) (x) dx ; 

& | (/>)yr-i = Jr(x)(^- 1 (/ ) )(x)dx. 
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It is worthwhile to point out, however, that, in absence of additional con- 
straints, neither JF -1 nor T is necessarily positive definite so that, for in- 
stance, the equation T | 0) = might have a non-trivial solution different 
from | 0) = 0. At this level, (0 | <j))p-i does not define a true norm but 
merely a pseudo-norm [17, 20]. 

It is evident that solving Eq. (2) amounts to analyzing the compatibility 
between C and H; in other words, C and H are to be found contextually. Once 
Eq. (2) is formally solved, one can investigate its supersymmetrization [21]. 
By this we mean the construction of super-charges 

«-(!!*), «-(£S) («> 

with anti-commutator 

k={q,q} = (^f (?) 

and a polynomial formulae 

n n 

TF = ]T a k H\ FT = £ al(H*) h (8) 

fc=0 fc=0 

with the final goal to elucidate the conditions leading to such a type of the 
closure of the algebra 



1.2 Plan of the paper 

In Section 2 we elaborate a particular solution to our problem inspired by 
the specific second-order supersymmetry (SSUSY) results of ref. [11]. Our 
solution of Eq. (8) will have the form 

FF = H*-j, FT = h\-°y (9) 

where h\ is naturally related to h 2 by Hermitian conjugation, 

h = h\ , (10) 

if c 2 is real. Our explicit solution to the problem is rendered possible by a 
SSUSY inspired gluing constraint [21, 22]. We show that 

Th 2 = h x T, (11) 

which, because of Eq. (10), is now equivalent to Eq. (2). This amounts to 

CVTh 1 = h 1 CVT . (12) 
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Explicit analytic examples of VT— asymmetric models are expressed in terms 
of circular or hyperbolic functions. 

In Section 3 we perform a detailed investigation of eq. (2) for a charge 
operator which is of the second order in derivatives, 

and where G (x) and D (x) are complex functions of the real coordinate x: 

G(x) = G R (x)+iG I (x), 
D(x) = D R (x)+iD I (x). 

We further derive the polynomial algebra of Eq. (8). In order to show 
explicitly that our formalism allows to generate VT— asymmetric models 
with real spectrum, we discuss in Section 4 a particular polynomial oscillator 
model. 

In Section 5 we generalize the postulate (13) and derive the general form 
of the charge operator C of any finite order in the derivative d/dx such that 
T = CV is Hermitian. At the very end, in section 6 we give some perspectives 
on the impact of our results on a variety of fields where the use of similar T 
might play significant role. 



2 SUSY gluing constraint 

Starting with a second-order C of the form (13) we have to guarantee, first of 
all, the Hermiticity of T = CV and = VC^ 1 . It is easy to show (see also 
section 5 below for an exhaustive discussion of these important conditions for 
polynomial charges) that the latter Hermiticity condition forces us to impose 
the necessary and sufficient requirements 

D R (x) = D R (-x) + ^G R (x), Dj(x) = -Dj(-x) + ^Gj(x) 

where G R (x) = G R (—x) is even while Gj (x) = —Gi (—x) must be odd. 



2.1 Factorization 

In the subsequent step of our considerations we factorize our second-order 
charge operator C as follows, 

C = qiq 2 , qi = ^- + U(x), q 2 = ^- + W(x), (14) 
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where 



U(x) + W(x) = G(x) , j-W(x) + U(x)W(x) = D(x) . (15) 

In order to simplify the problem at the start we impose the following "gluing" 
constraint on q\ and q 2 , 



= (q\)*qi + c, 



(16) 



where c is a complex number. By inserting Eqs. (14) into Eq. (16), we 
obtain 



i- + w] ( — i~+w 



dx 



dx 



dx ^ j \ dx j ' 



whence 



^-W(x) + W 2 (x) = —^U(x) + U 2 (x) + c. 
dx dx 



We find the following representation for TT* and T*T (Eq. (3)) 

= T (f)* = ( qi q 2 V) ■ (Vqtqiy = q x q 2 {V) 2 ((£)* (q§ 
which, taking Eq. (16) into account, becomes 



qi 



(qi)*qi + l + l (q\ 



2 2 

/ t\* c c 



qi 



c c 
2 ~ 2 



Correspondingly, 

,t 



Defining the Hamiltonian operators 



hi = qi(q\)* + 



^ ( t\* ^ c c 



'A 

' d 2 



U 



dx 



U 



c 
2 



(17) 
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and 



hi = T(qlYq 2 V-- 



d 2 d 



W(-x) + W 2 {-x) 



dx 2 dx 2 ' 

equation (7) provides the following representation for K 

k = h 2 -^, n = ( hl ^ 



4 ' \ h 
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Comparing with Section 3.3 below, where TT* = H 2 + aH + 7, and setting 
hi = H, we get in the present case a = and, correspondingly, Vq = 0, 
according to Eq. (53) below, as well as 7 = — c 2 /4. This shows explicitly 
how the present model can be derived from the general results of section 3. 



2.2 Hamiltonians 

Remembering the first of Eqs. (15), Eq. (17) becomes 

^-G(x) + W 2 (x) - (G(x) - W(x)) 2 = c, 

i.e., 

4~G(x) - G 2 (x) + 2G(x)W(x) = c, 
dx 

or 

G(x)W(x) = \ \G 2 {x) - -^G(x) + c ] . (18) 
We immediately deduce that 

G 2 (x)-±G(x) + c 
W x) = — — — — — 

, ^ G 2 (x) + ^-G(x) - c , N 

U(x) = G(x) - W(x) = V > ^ > . (19) 



Thus 



with 



v ™ v (G'(x)) 2 G"(x) G 2 (x) c 2 
VW = GW "W + 5GW + 4 + 4GW- (2 °» 
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From Eq. (18) we also get that at the zeros x of G, we must have 

— 

dx x=x 

which is a constraint on G, too. In fact, the method would fail if G had 
several zeros with non-identical values of the first derivative at each of them. 

An important comment must be made here since even if a function does 
not vanish on the real axis, one can investigate its zeros in the complex x 
plane. For instance, if 

1 + i sinh(ax) 



G(x) = G (x) = z(x) = 



a E R, 



(21) 



it is immediate to check that z(x n ) = at x n = — i (2n + 3/2)n/a, n = 0, 
±1, .... This would mean that dG (x n )/dx = (ia / 2) cosh(ax n ) = 0, thus 
implying that we must put c = in this case. 

In the similar spirit, we may consider the whole class of functions which 
depend on x only via z(x) of Eq. (21) in an arbitrary nonlinear manner, 
G m (x) = G(z(x)), since, as a function of x, z is PT-symmetric, and any real 
function of z is PT-symmetric, too, and is an acceptable candidate for G. 

It becomes convenient to change variables and express the Hamiltonian, 
H = —d 2 /dx 2 + V(x), with V(x) given by formula (20), as a function of z, 
by observing that 

dz d f— d 

2 



d_ 

dx 
dx 2 




2/1 

-or - 



d 2/1 
a z(l 

dz 



and 



V(z) 



d 




-a: 



z(l - z)^-G + a 2 

dz 



, z(l-z) d 2 CP e_ 
2G dz 2 + 4 + AG 2 



a 



] dz 2 ' 



,l-2z d 

AG Tz 



G{22) 



2.3 Consistency 

We prove now an important constraint on the complex number c = Cr + icj. 
From the second of Eqs. (15), we have 

d 



D(x) 



dx 



W(x) + U(x)W(x) 



2G 2 (x) 
1 



AG 2 



f 2G^-G - ^g) G-^-g(g 2 -^-G + c 
, dx dx 2 I dx \ dx , 



\dx I dx 
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or 



D{x) = 



D*(-x) 



AG 2 
1 

4G 2 



d 



d 



2G —G - 2G—G + —G + G 



dx 



dx 2 



dx 



4 c 2 



-2G 2 ^-G-2G^G+(^-g) +G 4 
ax ax 2 \dx ) 



(23) 



(24) 



where the functions on the right-hand-sides of Eqs. (23) and (24) are all 
computed at x. 

In deriving Eq. (24), use has been made of the fact that G and d 2 G/dx 2 
are PT-symmetric, while dG/dx is PT-antisymmetric, i.e., (dG / dx(—x))* = 
—dG/dx(x). Eq. (23) is obviously consistent with the general form of D as 
a function of G given by Eqs. (48), (50), with c 2 /4 = -I - D(x )G 2 (x ). 
Subtracting Eq. (24) from Eq. (23) side by side, we obtain 



D(x)-D*(-x) 



c 2 )* - c 2 



±G+^ 
dx 4G 2 



(25) 



Combining Eq. (25) with Eq. (33), we obtain the important result 

( c 2 )* - c 2 = -> 3(c 2 ) = -> c rCi = 0. 
From Eq. (19) we easily obtain 



U(x) = W*(-x) - 



CR 



G(xY 



whence 



and 



Thus 



{U*{x)f = W 2 {-x) + 



4 



(G*(x))< 



2c 



R- 



W(-x) 
G*(x) 



(26) 
(27) 



h\ — h 2 



^U(x) + U 2 (x)) + -^W(-x)-W 2 {-x)+c R 
dx I dx 



= c R 



(G 2 (x))* \dx 
Using Eq. (19) to replace W(—x), we obtain 
1 ( d 

a : - , s ... 

(G 



—G*(x) +c R \- 2 Y, / + 1 



G*{x) 



h\-h 2 = c R 



cr 



(G 2 )* \dx 
cr-c 



crCi . G 2 
% = -ic R c r . 



Therefore 



{G 2 Y U (G 2 )* U ^\G\ 4 ' 
h{ = h 2 CrCi = . 



(28) 
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2.4 Periodic potential 

Let us now give an example which generalizes PT-symmetric periodic po- 
tentials [8, 23]: 

G(x) = e iax + r , aGR, r G R, r ^ ±1. 

In this case we have, for all i£R, 

U(x) = l - (e iax + r) + l - [iae iax - c) / (e iax + r) 

= i (e im; + r) - ^ (iae im; - c) / (e im; + r) . 

Since G never vanishes, we do not have any constraint on the value of c, 
in addition to the one which requires that c be either real, or imaginary. 
The spectral analysis of the corresponding Schrodinger operators h\ and h,2 
with periodic potentials can be performed as a generalization to the non- 
"PT-symmetric case of the investigation done by in Ref. [24]. 

We now examine the invertibility of C and the boundedness of C _1 . First 
notice that C can be written in the following form 



where 



C = dC 2 , Cx=Cu + r -, C 2 = C v + r - (29) 

C u = -f + U 1 , U 1 = U- r - 
dx 2 

Cw = ~tL + Wu Wi = w- r -. 

We will discuss the invertibility of each factor in (29) separately. As for C\ 
we first observe that the numerical range {z = (Cuip,ip) '■ ip £ °f Cu 

is contained in the strip {z : \Rez\ < a} where a = max xe m\Ui(x)\. Hence, 
if |r| > 2a then — r/2 is in the resolvent set of Cu [25] and, therefore, C\ is 
invertible with bounded inverse on L 2 (R). A similar argument holds for Cw 
Thus, for sufficiently large values of |r|, operator C is invertible and C~ l is 
bounded on L 2 (R). 



3 Second-order charge operator C 
3.1 Re-construction of the potential 

We already noticed that in the second-order charge operator (13), the nota- 
tion of Section 5 below implies that we have the correspondences 72(2^) = 1, 
71 (x) = G(x) and jo(x) = D(x), so that the Hermiticity constraints on 
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the real and imaginary parts of jg(x), Eqs. (70) and (71), with u = 2 and 
£ = 0, 1, immediately give 

G R (x)-G R (-x) = 0; Gi (x) + d (-x) = 0; (30) 
D R (x)-D R (-x) = ^G R (x); D I (x) + D I (-x) = ^G I (x) .(31) 

As a consequence of Eq. (30), G is PT-symmetric 

G(x)=G*(-x), (32) 

while Eq. (31) yields 

D(x) - D*(-x) = 4-G(x) . (33) 
ax 

We assume that and H satisfy the intertwining condition (2) and that H 
depends on a local complex potential, V (x) : 

H = -£ + VW, (34) 

with V (x) = V R (x) + iVi (x). In turn, V R (x) and Vi (x) are conveniently 
decomposed into their even and odd parts: 

V R (x) = Vi(x) + Vg(x), 
Vj{x) = V I E (x) + V I °(x), 

with Vi (x) = Vi (-x) and Vg (x) = -Vg(-x),(K = R, I). We write now 
condition (2) explicitly and obtain three non-trivial equations by imposing 
that the coefficients of (d/dx) 2 , d/dx and (d/dx)° vanish, 

-2 (vg + iVi) + 2-^- (G R + iGj) = , (35) 



(36) 



dx 

2£ {Vi + iVg) - 2 £ (Vg + iVi) + 2± (D R + iDj) 
+£ (G R + id) - 2 (vg + iVi) (G R + id) = , 

(vi + ivp) - £ (vg + ivi) + £(d r + iDl ) 

+ (G R + id) ± (Vi (x) + iV? (x)) - (G R + id) £ (Vg + iVi) (37) 
-2{D R + iDi) (vg + tVi)=0 

while the coefficients of (d/dx) 4 and (d/dx) 3 are identically zero. 



dx 2 
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3.2 Integrability 

The first of the above equations (35) yields 

tf-E** v - E = l G - (38) 

The second equation (36)yields 

_ ^ _ ^ _ V o Gl + A_ D , + 1 £ G , = o, 

and is easily integrated for the other two components of the potential, V R (x) 
and Vf (x), as functions of G R (x), Gj (x), Dr (x) and Dj (x), by replacing 
V R (x) and Vf (x) with their expressions (38): 

Vi (x) = ~G R (x) + \ (Gr (*)) 2 - \ [Gj {x)f - D R (x) + V (39) 

V? (x) = l -^Gj (x) + G R (x) G: (x) - Dj (x) . (40) 

Here, Vq is a real integration constant. The corresponding integration con- 
stant in the equation for Vf(x) must be zero, because the function is odd. 
Both equations can be recombined as 

V(x) = ~G(x) + \g\x) - D(x) + V . (41) 

Finally, the third equation (37) allows us to express the Gj (x)'s, (J = R,I), 
as functions of the Dk (x)'s, (K = R,I), or, more conveniently, viceversa. 

-\j^ GR+q ih GR+ {t x Gr ) + K- g "- 2Dr )tx Gr 

Gi d 2 Gi ( d \ 2 . d d d 



, -Gj + 2 (Dj - GjG R ) —Gj - Gr—Dr + Gj—Dj 

2 dx z \dx J dx dx dx 

= (42) 

~Ydx^ + 2 Gr i^ Gi + y~ Gi +Gr ~ 2Dr ) Tx Gl + ~2 Gi 1x^ Gr 

+ 2^-Gr^-G: + 2(G R G I - Dj)^-G R - Gj^Dr - G R ^-Dj 
dx dx dx dx dx 

= 0. 

Eqs. (42) can be recombined in the following first-order linear equation 
expressing the unknown function D(x) in terms of the known function G(x) 
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and its derivatives up to third order 

Eq. (43) is easily solved by direct integration. Let us define the auxiliary 
functions 

g(x) = 2±G, (44) 

W)h (x) s wv (46) 

Eq. (46) is promptly integrated by use of definition (44) to 

pW . ex p(2/> GM )=IS), <47) 

where Xq is an initial point where G is different from zero. It is now easy to 
check that the general solution to Eq. (43) can be written in the form 

p(x)D(x) = [ X dx' ^/^ 1 +p(x )D(x ) , 

or 

The integral on the right-hand side of Eq. (48) is computed by elementary 
methods in the form 

G\x) G 2 (x)G'(x) G(x)G"(x) {G\x)f 

VbJu kj^Ju j j y^u ) — 

ixo 

with 



f dx'G(x')f(x>) = _ ^ W +^fL + I , (49) 



G 4 (x ) G 2 (x )G'(x ) G(x )G"(x ) (G'(x )) 2 
h = 4 2 + 2 — ' (50) 

where G' = dG/dx, and so on, thus providing the most general expression of 
D as a function of G and of its derivatives. 
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3.3 SSUSY algebra 

Assuming a charge operator, C(x), of the form (13), we now verify that the 
operator 

F{x)F*{x) = C{x)VC*(x)V = C{x)C*(-x) 
can be written as a particular case of formula (8) 

F(x)r*(x) =H 2 + aH + 1) 

where a and 7 are constants to be determined and H is Hamiltonian (34) 
with V given in (41). In fact, we have 

where use has been made of relations (30), (31) stemming from Hermiticity 
of C(x). After some algebra, the right-hand side of the above expression is 
brought to the form 

C(x)C*(-x) = ^ + (2D(x)-G 2 (x)-3G'(x))^- 2 (51) 

+ (2D'(x) - 3G"(x) - 2G(x)G'(x)) ^ 

+D"(x) - G"\x) + G(x)D'(x) 
-G(x)G"{x) + D 2 (x) - D(x)G'(x) , 

and is to be compared with 

H 2 + cH + ', = (-JL + V (x?J + a (-^ + v (z))+^ (52) 

-V"(x) + aV(x) + 7 , 

where V(x) may be expressed as a function of D(x) and G(x) according to 
Eq. (41). Direct comparison of the right-hand sides of the above formulae 
allows us to determine the a constant as 

a = -2V . (53) 

The value of 7 expresses the compatibility between C and the polynomial 
algebra through the equation 

V 2 (x) - V'\x)+aV(x)+ 1 = 
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= D"(x) - G"'{x) + G(x)D'(x) - G(x)G"(x) + D 2 (x) - D(x)G'(x). 

Here, we insert the expressions of V(x) and V"(x) in terms of G(x), D(x) 
and of their derivatives obtained from formula (41), and making use of Eq. 
(43), as well as of its general solution (48), (49), we obtain the final result 

1 = V 2 + I + D(x )G 2 (x ), (54) 

where I is defined in Eq. (50). This makes it possible to interpret 7 as 
a kind of integration constant. Thus, CVT invariance leads to the SSUSY 
polynomial algebra, Eqs. (7), (8). 



4 Polynomial oscillators 

The simplest factorization of C reads 



so that, correspondingly, 



D{x)= <m +G _M m 



In this case, Eq. (43) yields G"'(x) = 0, i.e., 

G(x) = ax 2 + ibx + c (57) 

where a, b and c are real numbers, owing to the fact that G(x) is VT- 
symmetric. From Eq. (41) we obtain: 

V(x) = l -G 2 (x) + G\x) + V 

11 11 c 2 

= -a 2 x 4 (b 2 - 2ac)x 2 + -iabx 3 + -xiibc + 4a) + ib H h V . (58) 

4 4 2 2 4 

If we make the additional assumption c = 0, for the sake of simplicity, the 
polynomial algebra provides the constraint 

on 7 [Eq. (54)]. 
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4.1 The problem of invertibility 

We will now make the spectral analysis for H and study the invertibility of 
T in the case c = 0. Then 

111 

V(x) = -aV - -b 2 x 2 + -iabx 3 + lax + ib + V . (59) 
4 4 2 

Setting /j 2 = ^ and v 2 = ^, we obtain an expression for the Schrodinger 
operator H of the same type as that presented in Eqs. (22), (23) of Ref. [19], 
namely 

d 2 

H = — — - + /iV - v 2 x 2 + 2ifivx 3 + 4/ix + 2iv + V (60) 
dx 2 

and D{H) = H 2 (R) n £>(:r 4 ), V/i, i/ G R, fi ^ 0. As in Ref. [19], H has 
discrete spectrum, i.e., the spectrum consists of a sequence of isolated eigen- 
values with finite multiplicity. 

In order to prove the reality of the spectrum of H, we first notice that H 
can be rewritten as 

d 2 

H = — — + x 2 (fix + iu) 2 + 4fj,x + 2iu + V . (61) 
dx 2 

Let us now perform the complex translation x — > x — Then if = S~ 1 H\S 
where Sip(x) = ip(x — on a dense set of functions ^ G -^ 2 (R) and 

d 2 



2ti> 



^ = 



/ w \ ( w ■ \ t 
+ \ x — — J ^//x — — + + 4/ix — liv + 2?z/ + V 



Hence H has the same spectrum of Hi. In turn Hi is selfadjoint on D(Hi) = 
D{H) = H 2 (R)r\D(x 4 ), thus it has real spectrum for all /i, v, V G R, fx ^ 0. 

We may stress that Hamiltonian (60) is not VT— invariant but has still 
a real spectrum because it is related by explicit similarity to the standard 
self-adjoint anharmonic oscillator. In our opinion this is an exceptional ex- 
ample since in general the proof of the reality of the spectra of non-Hermitian 
Hamiltonians cannot proceed in such a straightforward manner and, generi- 
cally, the necessary maps are non-local [26]. Moreover, by our construction, 
the reality of the spectrum is robust insofar as its CVT— symmetry cannot 
be spontaneously broken. In this sense, our example (60) may be perceived 
as a VT — asymmetric parallel to the VT— symmetric quartic oscillator of 
Buslaev and Grecchi [27]. 



• \ 2 / • \ 2 



16 



4.2 The problem of boundedness 



Let us now turn to the operator T = CV. In order to prove the invertibility 
of T and the boundedness of T~ x on L 2 (Tl) it is enough to demonstrate the 
same facts for C. Factorization (55) implies that it will suffice to prove that 
Ci = (J^ + f ) is invertible and that Cf 1 is bounded on L 2 (R) if G is given 
by (57). Indeed, we have 

d 1 i 
Ci = — + -ax 2 + -6x (63) 
dx 2 2 

and we now proceed as in Ref. [19]. More precisely 

^ d a ( ib\ 2 b 2 

C ^T. + 2{ X + Ta) + Ta (64) 

is similar to 

^ d a o b 2 

C * = T X + 2 X+ Za (65) 
via the complex translation x — > x — Hence Ci has the same spectrum 
as Ci- In turn C2 is unitarily equivalent, via the Fourier transformation, to 

^ ad 2 . b 2 

ft = -2^ + ,I+ 85- (66) 
Therefore C\ has the same spectrum as C3. Finally, we perform the unitary 
dilation (Uip)(x) = (a/2) 1 / 6 tp[(a/2) 1 / 3 x] and obtain that Ci has the same 
spectrum as 



,'n\ 1 / 3 

c, = 1 ■(-,(' ; : ' 



rf 2 . /a\-V3&2 

+ IX + - 



rfa; 2 V2/ 8a 



(67) 



Now, since the Schrodinger operator — -py + ix has an empty spectrum (see 
Ref. [28]), so does Ci. In particular z = belongs to resolvent set of C 1 , so 
that Ci is invertible and its inverse is bounded and defined on the whole of 
L 2 (R). 



5 Towards operators C of any finite order 

We shall postulate that the charge-operator component C of the pseudo- 
metric CV, where V denotes parity, is a polynomial of any finite degree 
lo — 0, 1, . . . in the momentum operator p, 

w d k 

C = E lk{x) ~~rr-> lk{x) = >y?(x) +iil{x) . (68) 
k=o ax 

The functions 7j^(x) and 7j£(x) are both assumed real, and our main task 
here is just to guarantee, at any integer oj, that the operator candidate for 
the metric CV is Hermitian. 
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5.1 The metric CV in differential form 

From 



k 



d {k-t) 



£ + m 



Rim)/ , - /(;//; 



df_ 

dx e 



(69) 



I(m), 



c f = E(-i) fc E 

k=0 1=0 
uj ( ui-£ 

= EHf E (-ir 

£=0 lm=0 

where the superscripts (m) at the functions 7 R and 7 7 indicate their m— tuple 
differentiation, one obtains that the Hermiticity condition CV = VC^ is equiv- 
alent to the (uj + 1)— plet of relations 



dx e ' 



uj-e 



VliV = l?(-x)+i-Yt(-x)= E I" 1 )' 



+ m 



m=0 



-R(m) / \ l ym ) / \ 

7f+m W — ? 7^+m \ x ) 



I(m) . 



with a trivial decoupling into its real and imaginary parts 



m 



if (-x) ~ l?{+x) = E ("l) m [ * Y l ) 



(70) 



and 



7/(-^) Wti+x) = - E ("l) m ( ) , (71) 

respectively, with £ = uj — k = 0,1, ... ,u. 

5.2 Functional freedom in complex coefficients 7fc(x) 

At the first few k — 0, 1, . . . the above Hermiticity constraints degenerate to 
the comparatively elementary relations, 

7*0*0 - 7«?(-z) =0, fc = 0, 
7*-i(*) - 7*-i(-*) = ( \ ) 7* (1) (*) , * = 1, 



a; 



7*- 2 (*) - 7*- 2 (-*) = ( W 1 J tH(*) v 2 
etc, or, in parallel, 



^ (2 \x), k = 2, 



lL(x)+li(-x)=0, k = 0, 

iL 2 (x) + ii- 2 (-x) -- 



V 2 
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etc. This means that the symmetric parts Hg(x) = Hg(—x) of all ^(x) are ar- 
bitrary functions while, in parallel, the antisymmetric parts hg(x) = —hg{—x) 
of all 7/ (x) are also arbitrary. We may conjecture that the remaining compo- 
nents Ri(x) = jg'ix) — He(x) = —Re(—x) and re(x) — jg(x) — he(x) = rg(— x) 
obey the rules 

Ru, = 0, R u -i(x) = < ^Hl 1 \x), R^ 2 (x) = ^±H i Jl 1 (x), ... (72) 
while 

r. = 0, r w _x(x) = ^hS\x), r^ 2 {x) = ^h^x), .... (73) 
and are fully determined by the respective recurrent relations (70) and (71). 



5.3 Proof 

We see that both the sequences R u ^ k {. x ) and r u _ k (x) have precisely the same 
structure so that just the sequence of R^-ki^x) may be considered without any 
loss of generality. Its elements should be evaluated in the recurrent manner 
with respect to the growing k. The appropriate Ansatze may be written in 
the finite-series form where, formally, H u+ i(x) = H u+2 {x) = ••• = and 
K+i(x) = K +2 (x) = • • • = 0, 

lS-k(x) = H„_ k {x) +Ec m > (74) 



rn=l 



(u-k)\ 



ll-ki*) = K-k(x) +Y. C m — ( k+ J^' h ( ^} k+m {x) . (75) 

m=\ \W - K)\ 

With an auxiliary cq — 1 these Ansatze describe all the u— dependence of 
our functions 7 = 7^ + i^ 1 in closed form. 

As already stated above, the first term and the subsequent sum are of an 
opposite parity in both these formulae since c 2n = at all n = 1, 2, . . .. This 
observation is easily proved since after the insertion of the latter two Ansatze, 
the complicated recurrences (70) are replaced by their simplified version 

2c 1 = -, 2c 2 = ---, 

% .e., 

^ = |(-D—' JJ^f. (™) 

It is worthwhile to point out that the c k coefficients with odd k can be written 
in terms of Bernoulli numbers (see, e.g., Ref. [29]) 

2(2 2n — 1) 

C2n-i= V (2yi) , } B 2 n (n>0). (77) 
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The key idea of an explicit solution of these recurrences is that the generating 
function f(x) = c k^ k of the coefficients c m must satisfy the functional 
equation f(x) — 2 = —f(x)/e x which is, in its turn, easily solvable. In this 
way we arrive at the solution of recurrences (76) in the following compact 
form, 

2 x 

fix) = c + cix + c 2 x 2 + . . . = r = 1 + tanh - = (78) 

1 + exp(— x) 2 

" 2(2»»-l) ,„_, 

Obviously, all the possible parity-violating terms in the right-hand side of 
our Hermiticity conditions (70) vanish. This makes the form of our polyno- 
mial charge C extremely flexible and confirms the consistency of its present 
construction. 

We may conclude that the requirement of Hermiticity of the metric CV 
defines all the antisymmetric components Re(x) and their spatially symmetric 
partners rg(x). It does not impose any additional constraint either upon the 
real and spatially symmetric coefficient functions H e (x) or upon their purely 
imaginary spatially antisymmetric partners hi{x). 



6 Outlook 

We now sketch some possible applications of our methods to a variety of 
problems where quasi-Hermitian or pseudo-Hermitian operators are involved. 

In the context of the Klein-Gordon description of the free motion of a 
spinless particle in the "usual" Hilbert space H the relativistic evolution is 
generated by the Feshbach-Villars [30] "Hamiltonian" H^ FV ) which proves 
non-Hermitian, 

\m) = e- iH ^-*o) mto)h H{FV) = _ l _ ^ 1 - A ^ (79) 

(cf., e.g., p. 341 in ref. [31]). One should notice that this model works 
with the differential pseudo-Hermitian operator with structure which strongly 
resembles the usual Schrodinger operators in the simplest non-trivial, two- 
dimensional coupled-channel case. Thus, we may expect that the methods 
described in our previous study might find an immediate extension to the 
similar problems. 

The idea may also find applications in a broader context, say, of the boson 
mappings in nuclear physics which were comprehensively discussed in the pa- 
per [13]. It is shown there that a consistent quantum mechanical framework, 
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and in particular a viable variational calculation for non-hermitian Hamilto- 
nians, can indeed be constructed after the introduction of a non-trivial met- 
ric. In the context of Holstein-Primakoff type mappings this freedom defines 
the link with so-called Dyson-Maleev type mappings (see [32]). In practical 
computations, a puzzling non-Hermiticity of observables proved more then 
compensated by the advantages, as has been amply demonstrated in appli- 
cations of generalized Dyson-Maleev mappings (see [33] and references cited 
therein) . 

All the technical conditions imposed upon the "true physical metric" G 
in review [13] are important, especially if one tries to work within a truly 
infinite-dimensional Hilbert space. This has been emphasized by Kretschmer 
and Szymanowski [14] who showed that the use of the toy metric operators 
might require a careful scrutiny because these operators remain unbounded. 
In this context, ref. [19] as well as our present paper demonstrated per- 
suasively that a switch to the use of the differential operators C might be 
understood as an important new idea. 

All the similar observations must be perceived as individual steps of a 
systematic improvement of the mathematically correct understanding of the 
use of the differential operators in connections with many applications of 
the quasi-Hermitian observables which seems to range, at present, from the 
elementary descriptions of the localization transitions in solid state physics 
[34] up to many ambitious VT— symmetric models in quantum field theory 
[35]. ' ^ 

The experience gained during our study of the simple Schrodinger equa- 
tions might equally well find applications on the very boundary of quantum 
mechanics (like, say, in cosmology [36]) or even in the domain of the classi- 
cal model-building (e.g., in the magneto-dynamics of fluids [37]) and in the 
various physical models of different origin characterized by the simple ma- 
trix structure of their description (see a number of their most elementary 
samples mentioned in the short and nice review [38]) where the eigenvalues 
coalesce or almost coalesce in the manner which contradicts the standard 
and robust finite-dimensional Hermitian-matrix mathematics. Of course, all 
these mathematical problems and not entirely standard physical situations 
may impose new and challenging tasks and motivate a deeper future analysis 
of the questions outlined in our present paper. 
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